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1. Introduction 

In the last few years Lattice QCD has become a precision tool, useful for the calculation of 
non-perturbative, gold-plated quantities in hadronic physics. This allows us to make stringent tests 
of QCD, and is crucial if we are to use Lattice QCD calculations as an ingredient in the search 
for physics beyond the Standard Model. It also makes it possible to determine accurately, from 
first principles, many parameters of the Standard Model (strong coupling constant, quark masses, 
elements of the CKM matrix.) 

Recently it has become possible to include charm quarks accurately in Lattice QCD calcula- 
tions, thanks to the introduction of the Highly Improved Staggered Quark (HISQ) action [jl|]. This 
is a fully relativistic discretization, and has a number of advantages, which we discuss later on, as 
well as some potential difficulties due to the much higher mass of c quarks as compared with u, d 
and s quarks. Our results for charm-light decay constants ph, as well as for the charm quark mass 
[Q] show that this approach can produce results with a similar accuracy to the previous staggered 
light quark calculations [Q]. 

Having seen that it's possible to use this approach with charm quarks, the question immediately 
comes to mind of whether the same is true for bottom quarks, with ensembles already existing or 
that can be realistically generated in the near future. In this exploratory study, we increase the mass 
of the HISQ quarks beyond the charm mass into the bottom mass region. 

We briefly discuss general issues with relativistic discretizations of heavy quarks in section 
|2| We then present in section || an update on our calculation of charm-light decay constants, and 
section § presents the status of a calculation of the ratio m c /m s , which when combined with our 
previous accurate value for m c will result in an improved m s determination. In section ||] we discuss 
preliminary results for masses and decay constants of bottom-bottom and bottom-strange systems, 
obtained through the relativistic discretization at masses beyond the charm quark mass and an 
extrapolation in the quark mass to get to the b. The same strategy is applied to the calculation of 
the b quark mass in section [6[ Section ^presents our conclusions and outlook for future work. 

2. Relativistic Discretization of Heavy Quarks 

The use of the ASQTAD improved staggered discretization in Lattice QCD has made it possi- 
ble, in recent years, to fully include the effect of u, d and s quarks on the QCD vacuum, resulting 
in realistic and accurate calculations of many gold-plated properties of hadrons [Q]. The ASQ- 
TAD discretization removes all tree-level a 2 errors present in the naive staggered action. Recently 
we introduced a new, further improved staggered action called HISQ (highly improved staggered 
quarks). The HISQ action goes beyond a 2 tree-level improvement by further smearing the gauge 
fields which enter in the staggered Dirac operator. We showed in ^ that this effectively cuts down 
the size of the 1-loop taste-changing interactions, which were typically larger than expected in the 
staggered formulation, as well as the ones liable to cause theoretical concerns if not under control 

A relativistic discretization has a number of advantages with respect to an effective heavy- 
quark theory such as NRQCD. If the discretized action has some symmetries, there are quantities 
which are free from renormalization. For example, for staggered quarks, pseudoscalar meson de- 
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cay constants do not renormalize because of PCAC. This eliminates one of the main sources of 
systematic errors from the calculation. Using the same formulation for the heavy and the light 
quarks is simpler, and allows us, for example, to calculate accurate ratios of quark masses. 

But using a relativistic action for heavy quarks brings forth new issues. Discretization errors 
will be generally proportional to powers of the quark mass, (am)", and the question is whether those 
can be kept under control for large masses. For staggered actions only even powers appear, and 
for an action like HISQ that removes all a 2 tree-level errors, the discretization errors start at (am) 4 , 
a s (am) 2 . These would still be too large for precision work; however, the fact that a system with 
charm (or heavier) quarks is essentially non-relativistic makes it possible, by a simple modification 
of the coefficient of the Naik term, to completely remove those errors at leading order in v/c, where 
v is the velocity of the quark inside the hadron Another essential ingredient is the existence of 
fine enough ASQTAD sea quark ensembles [0], going from a m 0.15 fm (very coarse), 0.125 fm 
(coarse), 0.09 fm (fine), 0.06 fm (superfine) and down to 0.045 fm (ultrafine). All these ingredients 
make it possible to use the HISQ action for calculations with charm quarks with almost the same 
level of precision as was possible for light quarks. The existence of ensembles at several lattice 
spacings is crucial to reliably extrapolate to the continuum limit. For bottom quarks, at present we 
must also extrapolate in the heavy quark mass: we do calculations at several values of the lattice 
spacing and the heavy quark mass (above the charm mass), and make a joint a — > 0, Mf, — ► Mb 
extrapolation. On the ultrafine ensemble such extrapolation is small. 



3. fo s decay constant 

Our prediction of the D and D s decay constants caused a lot of interest because experimental 
determinations soon after [[7], || [J gave a surprising picture: /o agreed very well with our result 



but fo, was more than 3 a higher [ ]10| , |1 Here a was dominated by the experimental error since 
our own error was so small. Because we use a relativistic discretization in which there is a PCAC 
relation, we can calculate and /b s from the partially conserved axial current, in exactly the same 
way as is done for f K and fg (for which we agree with the experimental values). Using three values 
of the lattice spacing, we obtained a 2% accurate result. We now have extended the calculation 
to two finer lattice spacings, the superfine a « 0.06 fm, and the ultrafine a « 0.045 fm. We show 
in Fig. [j] the preliminary results from this updated calculation, and a comparison of our data with 
current experimental determinations which have changed the picture considerably from that of last 
year. Now new CLEO determinations [12, [l3|] are closer to our result for fp s (especially for some 



channels) and a reanalysis of BaBar data [ ]14| ] is very close indeed. The discrepancy between our 
calculation and the current world average given by the Heavy Flavor Averaging Group is only about 
2a & \ 



4. Ratio of m c to m s 

Because we are using the same HISQ action for all quarks, light and heavy, we can obtain very 

'in fact the final discrepancy should be even a little smaller, because our final result will use the updated determi- 
nation of the parameter rl [p6p, which is used to set the scale and therefore the physical value of the lattice spacing, rl 
has decreased a little from its previous calculated value, and that in turn will increase slightly the value of /b s . 
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Figure 1: HPQCD lattice calculation and comparison with experiment. World average is from HFAG ^ 



accurate ratios of quark masses, where some systematic errors cancel. In particular, we have a very 
precise value for the ratio of lattice masses, m c /m s . We can now go to another scheme, like MS, 
and the Z factors that appear will cancel in the ratio; this allows us to combine this ratio with our 
precision calculation of the m c mass [j|] to obtain a determination of the strange quark mass with a 
« 1.5% total eiTor. 

In Fig. H we show preliminary results for m c jm s . The gray points in the plot are the lattice 
results for the ratio at various bare masses; the blue points are the results interpolated to the physical 
masses, determined by fixing the masses of the tj c and the 7] s to their physical values 2 ; the red point 
is the continuum and chirally extrapolated value. 

5. Masses and decay constants for bottom-bottom and bottom-strange systems 

In order to study bottom quarks with current lattices, we must resort to an extrapolation pro- 
cedure: we calculate at several values of the lattice spacing a and masses above the charm mass, 
Mh, and then we do a joint continuum and mass extrapolation, a — > 0, Mj, — > My. 

In Fig. |3| we show preliminary results for the difference between the masses of the heavy- 
strange and the heavy-heavy pseudoscalar mesons, AM/, S = M/ u —M^/l, for a range of masses 
from the charm to the bottom quark mass and for several lattice spacings. We can see that the results 
have a smooth dependence on both the heavy mass and the lattice spacing, and the extrapolated 
value to the bottom mass will be consistent with experiment. 

In Fig. H we show similar data for the decay constant of the heavy-strange pseudoscalar, f hs , 
which interpolates between fp s and fg s - The extrapolated value should be consistent, and already 



2 See Jl6| | for a discussion on what we mean by the "physical" value of tj s ; the main point is that it can be used to 
accurately fix the strange quark mass. 
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Figure 2: Ratio of charm to strange quark masses. 



competitive, with the previous NRQCD calculation [17]. 

6. Bottom quark mass 

Another thing which we can do with the relativistic formulation is a calculation of the bottom 
quark mass, with the same methods used for the charm quark mass calculation Let's recall that 
the basic idea is to take advantage of the high-order continuum QCD perturbation calculations per- 
formed by Chetyrkin, Kiihn, Steinhauser and Sturm, by combining them with high precision lattice 
data for moments of pseudoscalar and vector current correlators. For the case of a pseudoscalar 
current, the basic object is the correlator: 

G(t) = a 6 £ (am Qh ) 2 (0\j 5 (x,t)j 5 (0, 0) |0> , j 5 = Ws Wh (6-1) 

X 

which is finite and unrenormalized as a — > 0, due to PCAC. The quantities to compare with contin- 
uum PT are the moments, G n = £ r (?/a)"G(f), of such correlator, for which we have: 

" V ; (aM,^))' 1 - 4 

Here the g n are known to large orders in continuum PT for large Mf, and not too large n, and the 
left hand side is calculated from the lattice data, giving us the value of aMf,. In practice we use 
what are called reduced moments, in order to remove some of the discretization errors (see [§] for 
details.) 

In Fig. H we can see several reduced moments, for different values of the heavy-heavy pseu- 
doscalar mass and several values of the lattice spacing, as well as the continuum extrapolations. 
With this method we expect the final value on mj, to have a total error of 2 — 3%. 

7. Conclusions and Outlook 

It is feasible to use a highly improved relativistic action for charm quarks on fine enough 
lattices. This allows us to do very accurate calculations for systems with charm quarks. We obtain 
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Figure 3: AM hs = M hs - M hh /2 



a 2% accurate result for fo,. Due to changes in the experimental values, our calculation is now 
only about 2a away from experiment. We also can obtain a very accurate m c /m s mass ratio, which 
combined with our 1% m c determination will provide a sa 1.5% determination of m s . 

To get to the bottom quark mass, we need (at least for now) to extrapolate in the mass of the 
heavy quark. This can be done, as we show for AM^ = Mf, s — M^/2 and fg s , and it seems that the 
additional source of error that it introduces can be kept under control. The results are in both cases 
consistent with either experiment or previous NRQCD calculations. 

If we had even finer lattices (0.03 fm, say) we could reduce substantially that uncertainty. This 
may be a good strategy to reduce the errors on some B physics quantities, for example fg s . 
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